Photonic band gap via quantum coherence in vortex lattices of Bose gases 
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We investigate the optical response of an atomic Bose-Einstein condensate with a vortex lattice. 
We find that it is possible for the vortex lattice to act as a photonic crystal and create photonic band 
gaps, by enhancing the refractive index of the condensate via a quantum coherent scheme. If high 
enough index contrast between the vortex core and the atomic sample is achieved, a photonic band 
gap arises depending on the healing length and the lattice spacing. A wide range of experimentally 
accessible parameters are examined and band gaps in the visible region of the electromagnetic 
spectrum are found. We also show how directional band gaps can be used to directly measure the 
rotation frequency of the condensate. 

PACS numbers: 42.50.Gy, 42.70.Qs, 74.25.Qt 



In-situ optical probing of vortices in rotating atomic 
Bose-Einstein condensates (BECs) is a challenging task 
due to the small size of vortices requiring imaging resolu- 
tion close to a single wavelength. It is usually necessary 
to let the condensate expand ballistically before employ- 
ing imaging methods [lj, |2|, y, Li yf • In-situ imaging of 
the vortex lattice has been demonstrated, by imaging the 
cloud perpendicular to the axis of rotation and observ- 
ing patterns created by the overlapping lattice planes Q . 
Alternative schemes based upon ultra-slow light imaging 
have been proposed for non-destructive, in-situ dynami- 
cal imaging of topological and collective excitations0 as 
well as for probing phase evolution of the condensate 

In a simple treatment of light propagation in a vortex 
lattice of a BEC, one may ignore the spatial periodicity 
in the system as the vortex cores and the atomic sam- 
ple would have negligible difference in their refractive in- 
dices. The present paper however, aims to exploit this 
periodicity to generate photonic band gaps, similar to 
those in photonic crystals which are well-known for their 
light localization effects By modifying the dispersion 
relation of light, photonic crystals create many optical 
effects ranging from left-handed light propagation [lOj to 
slowing down or stopping light [llj . If the vortex lattice 
of the condensate may be utilized as a photonic crystal, 
these optical effects can be used to modify or probe the 
condensate, or create novel atom-optical devices. 

As a first example, presence of a photonic band gap 
may provide stability to the vortex lattice against ra- 
diative losses Such losses are especially problem- 
atic as a source of decoherence in the quantum in- 
formation applications, based upon generation, stor- 
age and transfer of nonclassical states in atom-optical 
systems |T3, fill ITEl fl6l \v\ . These operations typically 
involve excited states and would suffer from spontaneous 
emission. The idea of suppressing spontaneous emission 
and enhancing stability of a Bose condensed system us- 
ing photonc band gaps was first discussed for the case of 
optical lattices [T^. Unfortunately, due to the competing 



effects of weak excitations and the strong periodic po- 
tential which are both provided by the same laser beams 
forming the optical lattice, the suppression was found to 
be insignificant 18]. In the case of a vortex lattice, how- 
ever, the periodicity is a topological effect and the atoms 
can be excited by a completely different and indepen- 
dent physical system such as a pump laser. Hence, one 
would expect to achieve a more significant suppression 
in the case of vortex lattices, which could be helpful to 
improve performance of atom-optical devices. One may 
also exploit the photonic band gap as an alternative to 
the existing EIT schemes for slowing down or stopping 
light, as discussed in a number of recent papers pT| . 

In this letter, we discuss presence of photonic band 
gaps in a vortex lattice and show how these band gaps 
can be used to measure the rotation frequency of a vor- 
tex lattice. In vortex lattice experiments, the rotation 
rate of the vortex lattice is not directly measured, but 
rather inferred from the change of the cloud's aspect ra- 
tio upon rotation^], 0, 0, 0, U • It is assumed that, with 
a fully developed vortex lattice, the cloud rotates like a 
solid body; on the other hand coupling of the collective 
modes of the lattice with the cloud's collective oscilla- 
tions may result in a much more complicated behavior. 
An independent measurement of the rotation frequency 
of the lattice would clarify the dependence of the cloud 
aspect ratio on rotation as well as determine clearly how 
far one can probe into the lowest Landau level regime 
0- 

In order to have a significant photonic band gap, it is 
preferable to have a material consisting of periodically 
alternating sequence of dielectric domains with high con- 
trast in their dielectric constants. Dilute atomic gases 
exhibit strong differences from the dielectric constant of 
free space only at their resonance frequencies, where the 
optical fields would suffer high absorbance. To overcome 
this challenge, we consider utilization of a quantum co- 
herent index enhancement scheme which allows for many 
orders of magnitude improvement in the refractive in- 



dex of an atomic gas at frequencies where the absorption 
can be made zero[2(|. This effect is based upon the simi- 
lar principle of quantum coherence in electromagnetically 
induced transparency (EIT) [2l|,which is used for slow- 
light [22j and frozen light effects 23] . 

There are various schemes in which index enhancement 
accompanied by vanishing absorption is possible, such as, 
the initial coherence scheme, the Raman driven scheme, 
and the microwave driven schemes [3 H3- As an exam- 
ple, the upper-level microwave scheme is illustrated with 
the energy level diagram in the inset of Fig. ^ Quan- 
tum coherence is established through coupling of levels a 
and c via a strong microwave field, which allows for ab- 
sorption free propagation of a probe field dipole coupled 
to levels a and b. A pump mechanism from level b to c 
ensures an ultra high refractive index by maintaining a 
small fraction of atoms in level c. At the point of vanish- 
ing absorption, susceptibility (%) becomes a real number 
and can be estimated to be x ~ 3iVA 3 /407r 2 , where N 
is the atomic density, and A is the wavelength of a — b 
transition. In a dense BEC, the first correction to the sus- 
ceptibility is argued to be the local field correction^, 
which changes x to x/(l — x/3). Local field corrections 
are small below N ~ 10 21 m~ 3 at which \ increases to 
3.22 from 1.55 for 23 Na gas of A = 589nm. For 87 Rb 
gas of A = 794nm, N — 5.5 x 10 20 m -3 , we find x ~ 7 
including the local field correction. The dispersive be- 
havior of e near the frequency of index enhancement will 
be ignored. This seems reasonable as x is slowly varying 
about its peak value where the first derivative of x van- 
ishes and the first correction is of second order in probe 
detuning. 

As our atomic medium, we consider a vortex lattice 
of rotating BEC. It has been predicted d EHH El, 
and observed[H S H S H 13 that the vortex lattice is 
triangular. It has recently been shown that triangular 
lattice is robust under variations of the trap potential in 
two dimensions j3l|. The number of vortices can be as 
large as few hundred. Presently, we assume the lattice 
is large and ignore finite size effects. Lattice spacing a 
can be compared to the healing length £ = 1/ ^^SnNa^, 
where a sc is the scattering length. We consider a ~ 5 — 
10£. Our examinations indicate a rich variation of band 
structure depending on a. Below, we report photonic 
band gaps in the visible range. 

In the two-dimensional geometry, we consider straight 
line vortices with perfect cylindrical symmetry for sim- 
plicity. For vortices as effective cylinders of radius 
R = 2£, this corresponds to a filling ratio / = 
(27r/\/3)(i? 2 /a 2 ) ~ 14%. To study the effects of the vor- 
tex profile on the band structure, we also consider a sim- 
pler system of one dimensional vortex lattice. Such vor- 
tex row structures have theoretically been predicted|3ll 
|3l, but have not been observed so far. We use a vortex 
profile p(z) in the Pade approximation (3^. Our calcu- 
lations show that the form of vortex profile causes little 
change in the band structure. Assuming a dilute thermal 
gas background with uniform density, dielectric constant 
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FIG. 1: (Left) Position dependence of the dielectric function 
of the one dimensional vortex lattice within a spatial period 
where the background dielectric constant is assumed to be 
eo = L Inset shows the level scheme of the upper- level mi- 
crowave driven scheme. (Right) Band structure of a one di- 
mensional vortex lattice. 



within the vortex core region will be taken as eo = 1. 
Defining e\ = 1 + x, dielectric function of the atomic 
system can be expressed as e(z) = e a + (ei — e )p(z) = 
1 + xp( z )- We calculate the band structure of the vortex 
lattice of N = 6.6 x 10 20 m~ 3 23 Na atoms using plane 
wave expansion method with 11x11 grid of plane waves. 
Our results for a — 10£ are shown in FigQ There is 
a full gap at uja/2Trc = 0.35 with a gap-midgap ratio 
of 0.143. In a similar band structure for 87 Rb, with 
N = 5.5 x 10 20 m~ 3 , full gaps at visible midgap frequen- 
cies 520THz and 445THz are found when a — 10£ and 
a = 5£, respectively. 

Typical results for TE polarization are shown in FigEl 
for the two-dimensional lattice. Similar band structure 
occurs for TM polarization. In the TM case, the gap be- 
tween the first and second bands decreases and becomes 
negligibly small at the K point of the Brillouin zone. For 
the TE case, a gap appears at all points in the irreducible 
Brillouin zone including the M point. There are direc- 
tional pseudogaps but no complete band gap for both TM 
and TE polarizations. It should be emphasized that it is 
actually possible to observe full band gap, as it is typical 
for triangular lattice, for another choice of parameters. 
Assuming a >~ 4£ gives R/a <~ 0.5 for the radius of 
the vortex core R = 2£. A full gap for TE polarization 
is found for 87 Rb with a = 5£ at N = 5.5 x 10 20 m~ 3 . 
Full gaps for both polarizations may appear at higher 
index-contrasts. 
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FIG. 2: TE band structure of a triangular vortex lattice with 
a = 10£ oi N = 6.6 x 10 20 m~ 3 23 Na atoms. Different line 
patterns are used for visual clarity. 
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FIG. 3: (Color online) A vortex lattice with 31 vortices. Vor- 
tex spacing is taken to be 5 times the vortex core size which is 
equal to the light wavelength. Vortex cores have n = 1 while 
BEC has index of refraction n — 3. With these parameters 
infinite lattice results indicate a directional band gap, which 
can be observed by comparing this figure with the Fig0] Di- 
rection of the incident light is denoted by an arrow in both 
figures. 



After investigating an infinite vortex lattice and find- 
ing that band gaps develop, we consider whether the ef- 
fects discussed before can be observed in an experiment 
where the vortex lattice is necessarily finite. In a typ- 
ical experiment, there are three major differences from 
the idealized vortex lattice considered above. Here, we 
consider all three separately. 
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FIG. 4: (Color online) Same finite vortex lattice as in FigHTl 
but rotated by n/6. A directional band gap is clearly seen as 
there is no transmission. 



First, we assumed a strictly two dimensional vortex 
lattice, which disregards any bending or motion of the 
vortices along the rotation axis. However, the bending 
of the vortices have theoretically been predicted to be 
small except near the edges of the cloud 0, 13 ■ This 
prediction was later confirmed by viewing the vortex lat- 
tice from the side @. As the deviation of the position 
of the vortices are much less than the vortex spacing, we 
expect the bending, or Kelvin mode excitations of the 
vortices in the third dimension to be a very small effect. 

Second simplification in our calculation was that the 
index profile of the cloud remained uniform except for 
the appearance of vortices. In a trapped BEC with vor- 
tices overall density profile is maximum in the middle and 
would slowly decay towards the edges. Average density 
profile is predicted to be Gaussian or Thomas-Fermi 
like |35j | , in both cases a monotonous slowly varying func- 
tion. This slow variation of the index profile would create 
an extra source of scattering from the cloud, especially 
under coherent index enhancement. The presence of this 
extra scattering may complicate the optical band struc- 
ture effects if the variation of the density is not slow 
on the scale of vortex spacing. However, as the over- 
all density decays monotonously, a fast variation of the 
density on the scale of vortex spacing would mean there 
are only a small number of vortices in the cloud. So for 
a BEC containing enough number of vortices to create 
band structure effects, photonic band structure should 
be superimposed on the extra scattering effect which is 
expected to vary slowly in frequency. Thus, the photonic 
crystal effects should be easily observable. 

A third, and related to second, concern in an exper- 
iment would be the finiteness of the vortex lattice. We 
performed a numerical simulation of the propagation of 
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an electromagnetic wave in a finite lattice to quantify the 
number of vortices needed to observe band structure ef- 
fects. In Figs l3l4l we display our results for 31 vortices in 
two different orientations. The simulation was done with 
parameters where the infinite lattice calculation indicates 
a directional band gap, and this band gap is clearly seen 
by contrasting the two figures. In FigEl a Gaussian beam 
passes through the vortex lattice without much scatter- 
ing, while in FigQ] after a rotation of the vortex lattice by 
7r/6, there is no transmission. Our results indicate that 
around 20 vortices are enough to qualitatively observe 
the infinite lattice predictions for the band gaps. 

An important consequence of such a directional band 
gap is that it can be utilized to measure the rotation fre- 
quency of the BEC. In the vortex lattice experiments, 
the rotation frequency of the vortex lattice is measured 
indirectly by inferring it from the radius of the cloud 
[B S S U [j 0] However, if a directional band gap is 
established, probing laser beam will be chopped by the 
rotation of the vortex lattice and counting the chopping 
frequency will yield a direct measurement of the rota- 
tion frequency. Such a measurement can be carried out 



non-destructively, and continuously. As the rotation fre- 
quency of the cloud, which is around hundred Hertz in 
a typical experiment, is much smaller than the laser fre- 
quency, rotation will not affect the optical band structure 
discussed earlier. 

We have shown that photonic band gaps can be gen- 
erated in vortex lattices of an atomic Bose-Einstein con- 
densate by utilizing quantum coherence based ultra-high 
index enhancement effect. We considered a two dimen- 
sional triangular lattice and found the physical parameter 
regime for band gaps and other band structure effects to 
be observed. We argued that the effects found for an 
infinite lattice could be observed in experiments where 
tens of vortices are created in the condensate. We fi- 
nally suggested a method to use a directional band gap 
to measure the rotation frequency of the condensate non- 
destructively. 
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